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ABSTRACT 



Several studies have attempted to characterize morphological brain changes due to chronic pain. Although 
it has repeatedly been suggested that longstanding pain induces gray matter modifications, there is still 
some controversy surrounding the direction of the change (increase or decrease in gray matter) and the role 
of psychological and psychiatric comorbidities. In this study, we propose a novel, network-oriented, meta- 
analytic approach to characterize morphological changes in chronic pain. We used network decomposition 
to investigate whether different kinds of chronic pain are associated with a common or specific set of altered 
networks. Representational similarity techniques, network decomposition and model-based clustering were 
employed: i) to verify the presence of a core set of brain areas commonly modified by chronic pain; ii) to 
investigate the involvement of these areas in a large-scale network perspective; iii) to study the relationship 
between altered networks and; iv) to find out whether chronic pain targets clusters of areas. Our results 
showed that chronic pain causes both core and pathology-specific gray matter alterations in large-scale 
networks. Common alterations were observed in the prefrontal regions, in the anterior insula, cingulate 
cortex, basal ganglia, thalamus, periaqueductal gray, post- and pre-central gyri and inferior parietal lobule. 
We observed that the salience and attentional networks were targeted in a very similar way by different 
chronic pain pathologies. Conversely, alterations in the sensorimotor and attention circuits were differentially 
targeted by chronic pain pathologies. Moreover, model-based clustering revealed that chronic pain, in line 
with some neurodegenerative diseases, selectively targets some large-scale brain networks. Altogether these 
findings indicate that chronic pain can be better conceived and studied in a network perspective. 

© 2014 Published by Elsevier Inc. 
This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Chronic pain is defined by the International Association for the 
Study of Pain as a pain persisting over the healing phase of an injury 
(Loeser and Treede, 2008). Whether continuous or recurrent, chronic 
pain must be of sufficient duration and intensity to adversely affect a 
person's level of function, well-being and quality of life, chronic pain 
states are generally considered as either neuropathic or nociceptive, 
although this subdivision has been criticized and less clear-cut sub- 
divisions are proposed (Bennett et al., 2006; Treede et al., 2008). 

Converging evidence from animal and human studies indicates 
that chronic pain induces a dramatic anatomical and functional reor- 
ganization of brain structures and networks (Apkarian et al, 2009). 



Studies using voxel-based morphometry (VBM), a computational- 
based technique that measures focal differences in concentrations 
of brain tissue using a voxel-wise comparison between two groups 
of subjects (Ashburner and Friston, 2000), have reported conflicting 
results suggesting that chronic pain can either decrease (Apkarian 
et al, 2004; Buckalew et al, 2008; Burgmer et al, 2009; Davis et 
al, 2008; Draganski et al, 2006; Geha et al, 2008; Gerstner et al, 
2012; Gustin et al, 2011; Gwilym et al, 2010; Kuchinad et al, 2007; 
Rocca et al, 2006; Rodriguez-Raecke et al, 2009; Schmidt-Wilcke et 
al, 2006; Schmidt-Wilcke et al, 2010; Seminowicz et al, 2011; Tu et 
al, 2010; Unrath et al, 2007; Valet et al, 2009; Valfre et al, 2008; 
Vartiainen et al, 2009), or increase (Etgen et al, 2005; Garraux et 
al, 2004; Obermann et al, 2007; Obermann et al, 2009; Schmidt- 
Wilcke et al, 2006; Schmidt-Wilcke et al, 2007; Schweinhardt et al, 
2008; Seminowicz et al, 2011; Tu et al, 2010; Younger et al, 2010) 
gray matter density. Some authors have even reported a lack of any 
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change (Hsu et al, 2009; Rocca et al, 2006; Schmidt-Wilcke et al, 
2005; Schmidt-Wilcke et al, 2007) in gray matter volume or density. 
These discrepancies may be underpinned by the different category 
of chronic pain considered (e.g. low back pain or migraine), by the 
different underlying etiology (e.g. nociceptive vs neuropathic) or by 
other confounding factors such as, for example, concomitant psychi- 
atric comorbidities like depression or anxiety. Many fields of cognitive 
neuroscience are now moving towards a network approach to capture 
the dynamics of brain functioning (Bressler and Menon, 2010). New 
views posit that the elaboration of noxious stimuli and pain percep- 
tion depend largely on the activity of networks (Cauda et al., 2013; 
Mayhew et al., 2013). Recent findings by our and other groups have 
indicated that noxious stimuli trigger the response of different brain 
networks at different temporal delays from stimulus onset (Cauda et 
al., 2013; Mayhew et al., 2013). Network dynamics have also been 
proposed to underpin the process of chronification of pain (Baliki et 
al, 2012; Farmer et al, 2012; Hashmi et al, 2013). However, no study 
has so far characterized gray matter alterations in chronic pain using 
a network approach. In the present study, we performed a meta- 
analysis of VBM studies investigating structural changes in the gray 
matter of chronic pain patients using a novel network approach. We 
addressed three questions: i) is there a core set of brain areas altered 
in chronic pain?; ii) are different kinds of chronic pain associated with 
a common or specific set of altered large-scale brain networks?; and 
iii) when alterations of gray matter are reported together, do they 
reflect clusters of modified areas that share some relationship or are 
they casually found together? 

2. Materials and methods 

2. 1 . Selection of studies 

Systematic searches using 'chronic pain' and 'voxel-based 
morphometry' as keywords were conducted on the Medline, Sco- 
pus and Scirus databases. Additional searches were performed for 
the term 'voxel-based morphometry' associated with each disor- 
der indicated by the American Chronic Pain Association (ACPA) as 
part of the spectrum of chronic pain (http://www.theacpa.Org/7/ 
Conditions.aspx). Reference lists of the identified papers were ex- 
amined for studies not found with the database search. As a result 
we selected 80 papers. Of these, we excluded: a) non-original stud- 
ies (N = 17); b) studies that did not include gray matter location in 
Talairach/Tournoux or in Montreal Neurological Institute (MNI) co- 
ordinates (N = 6); c) studies in which the field of view was confined 
to a restricted region of the cortex (N = 13); d) studies in which the 
comparison between chronic pain patients and healthy subjects was 
lacking and e) studies with unspecified VBM analysis (N = 1 0). We also 
tried to identify any instances of multiple reports of single data-sets 
across articles, to ensure that only one report of a study contributed to 
the coordinates for the present meta-analysis (see PRISMA, 2009 Flow 
Diagram of article selection in Supplementary materials for further 
details). 

Following the recommendations provided by Rainville and Duncan 
(2006), we used tables to collect a detailed description of neuroimag- 
ing modalities, VBM techniques and data reporting the number of 
increased and decreased gray matter foci. We also checked the clin- 
ical condition of the experimental samples in terms of adherence to 
the diagnostic criteria for chronic pain, any comorbidity, medication 
and the overlap of control groups and experimental samples. 

Based on these criteria, 32 papers were included, with a total 
of 1509 subjects (739 patients and 707 controls). There is a lack of 
information concerning the gender of the sample for three studies 
(Schmidt-Wilke et al, 2005; Oberman et al, 2007 Study 1 ; Buckalew 
et al. 2008). Given all the others, the patient group comprised 196 
men and 553 women. The control group comprised 252 men and 565 
women (Table 1 ). 
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2.2. Anatomical likelihood estimation (ALE) 

We used the anatomical likelihood estimation (ALE) procedure to 
evaluate the presence of a common pattern of gray matter alterations 
in all types of chronic pain. ALE meta-analysis is a quantitative voxel- 
based meta-analysis method that can be used to estimate consistent 
activation (or areas of gray matter increases/reductions) on the basis 
of foci of interest across different imaging studies that have reported 
statistically significant peaks of activation (Laird et al, 2005b; Laird 
et al, 2009). During an ALE analysis, each activation focus is mod- 
eled as the center of a Gaussian probability distribution to generate 
a modeled activation (MA) map for each reported study. These 3D 
Gaussian distributions are consequently summed to generate a sta- 
tistical map that estimates the probability of activation for each voxel 
as determined by the entire set of studies. This map is then thresh- 
olded using a permutation test (Laird et al., 2005a; Lancaster et al, 
2000; Lancaster et al, 2007). 

We used a new ALE algorithm that estimates the spatial un- 
certainty of each focus taking into account the possible differences 
among studies related to sample size (Eickhoff et al, 2009; Eickhoff 
et al, 2012). This algorithm comprises a method to calculate above- 
chance clustering between experiments (i.e. random effects analysis, 
RFX), rather than between foci (fixed effects analysis, FFX) (Eickhoff 
et al, 2009; Eickhoff et al, 2012). In view of this, when we gathered 
data that were non-Talairach coordinates, we performed an accurate 
transformation using the most recent and unbiased method (Eickhoff 
et al, 2009). 

ALE maps were computed using a Java-based version of the ALE 
software named GingerALE (version 2.0.4) and customized Matlab 
routines at an FDR-corrected threshold of p < 0.05 and a minimum 
cluster size of K > 50 mm 3 . 



2.3. Jackknife analysis 

To rule out the possibility of some activations being driven by the 
involvement of a small subset of studies, we performed a jackknife 
analysis. The jackknife is a non-parametric method for estimating the 
sampling distribution of a statistic (Radua et al, 2011; Radua and 
Mataix-Cols, 2009). Given a sample data-set and a desired statistic 
(e.g. the mean), the jackknife works by computing the desired statistic 
with an element (or a group of elements) deleted. This is done for 
each element of the data-set. These statistics are used to generate an 
estimate of the sampling distribution. The results are presented as 
probability maps where a high probability means that a certain area 
is led by the majority of the experiments included. By contrast, a low 
probability means that a certain area is led by few experiments thus 
indicating the need for caution when interpreting this point. 



2.4. Network decomposition 

In a previous study, Biswal and colleagues used a large cohort of 
volunteers (1414) who underwent a resting-state fMRI scan to par- 
cellate the brain surface, and found that, with the use of data-driven 
methods, the resting brain can be clustered into 20 large-scale net- 
works, which are also identifiable when the brain is involved in an 
active task (Biswal et al, 2010). We applied the same panellation 
technique to the selected studies and examined how many voxels of 
altered gray matter density fell within each network. This allowed 
us to investigate whether chronic pain is associated with differential 
modifications of brain structures belonging to different brain net- 
works, such as attentional, thalamic and sensorimotor brain struc- 
tures. 
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Table 1 

Overview and demographics of the included studies. 

Sample Healthy control 

Range Range 





Year 


First author 


Subjects Age(Y) 


(Y) 


/ 


Diagnosis 


Duration (Y) 


Subjects Age(Y) (Y) 


/ 


— ni TT- 

Neuropathic 


201 1 


Gerstner 


9 


25.5 ± 2.5 




0/9 


Myofascial tem- 


2.5 ± 2.1 


9 


24.8 ± 1.4 


0/9 


pain 














poromandibular 
























disease 












2011 


Gustin 


21 


45.7 ± 2.9 


28-70 


4/17 


Trigeminal 


8.5 ± 2.1 


30 


53.6 ± 3.2 24-87 


6/24 
















neuropathy; 
























trigeminal 
























neuralgia 












2010 


Younger 


14 


38 ± 13.7 


23-61 


0/14 


Myofascial tem- 


4.4 ± 2.9 


15 


38 ± 13.7 23-61 


0/15 
















poromandibular 
























disease 












2010 


Schmidt- 


8 


52.2 ± 8.9 




1/7 


Persistent 




11 


51.3 ± 8.6 


2/9 






Wilcke 








idiopathic facial 






















pain 












2009 


Vartiainen 


8 


47 


41-53 


V7 


Labial or genital 


3-20 


28 


32 22-53 


10/18 














herpes 










2009 


wood 


30 


42.03 ± 8.43 




0/30 


Fibromyalgia 




20 


40.05 ± 10.01 


0/20 




2008 


Geha 


22 


40.7 ± 2.3 




3/19 


Complex regional 


3-13.5 


22 


40.5 ± 2.3 


3/19 
















pain syndrome 












2008 


Buckalew 


8 






4/4 


Chronic back pain 


3 months + 


8 


69.9 ± 3.9 






2007 


Schmidt- 


20 


53.6 ± 7.7 




1/19 


Fibromyalgia 




22 


50.7 ± 7.3 


2/20 






WllCKe 






















2007 


Kuchinad 


10 


52 




0/10 


Fibromyalgia 




10 


45 


0/10 




2006 


Draganski 


28 


41.9 ± 13.8 




20/8 


Phantom limb pain 


119.3 ± 140 m 


28 


41.4 ± 13.7 


20/8 




2006 


Schimdt- 


18 


50.4 ± 6.8 


34-59 


9/9 


Chronic back pain 


176 ± 87.60 m 


18 


49.9 ± 8.7 32-60 


9/9 






Wilcke 


















2004 


Apkarian 


26 






10/16 


Musculoskeletal; 


11.3 ± 11.4 


26 


43.3 


10/16 














radiculopathy 








Nociceptive 


2011 


Barnden 


25 


32 


19-46 


6/19 


Chronic fatigue 


7.4 ± 3.5 


25 


32.8 20-46 


6/19 


pain 














syndrome 












2011 


Pantano 


19 


53.2 ±11.2 




4/15 


Primary cervical 


12.7 ± 6.5 


28 


47.5 ± 15.6 


11/17 














dystonia 










2010 


Seminowicz 


56 


32.2 ± 12.3 




0/56 


Irritable bowel 


11.1 ± 7.7 


49 


31.1 ± 12.3 


0/49 














syndrome 










2010 


Gwilym 


16 


68 




8/8 


Primary 




16 


68 


8/8 
















osteoarthritis 












2010 


Tu 


32 


23.84 ± 2.99 




0/32 


Primary 


10.19 ± 3.25 


32 


23.81 ± 2.8 


0/32 
















dysmenorrhea 












2009 


Obermann 


7 


39.6 ± 15.1 


18-65 


0/7 


Chronic 


118.1 ± 51.4 


30 


35 ± 13.7 19-67 


13/17 
















posttraumatic 
























headache 












2009 


Rodriguez 


32 


66.8 ± 9 




13/19 


Primary 


7.35 


32 


63.9 ± 8.8 


13/19 
















osteoarthritis 












2009 


Valet 


14 


51.1 


28-68 


0/14 


Miscellanea 


9.8 ± 7.2 


25 


51.7 32-60 


11/14 




2008 


Kim 


20 


33.7 ±11.3 


15-53 


3/17 


Aura 


9.8 ± 6 


33 


33.8 ± 10.5 15-53 


4/29 




2008 


Schweinhardt 


14 


25.7 ± 5.1 


19-36 


0/14 


Vulvar vestibulitis 


5 ± 2.9 


14 


25.6 ± 6 


0/14 














syndrome 










2008 


Valfre 


11 


38.9 ± 6.4 




2/9 


Aura 


zU.b ± o.y 


27 


34.9 ± 8.6 


7/20 




2008 


Davis 


9 




30-58 


3/6 


Irritable bowel 




11 


24-50 


4/7 














syndrome 










2007 


Unrath 


63 


63.7 ±11.4 




18/45 


Idiopathic restless 


22.3 ± 7.1 


40 


63.4 ± 9.9 


1 1 /29 
















legs syndrome 












2007 


Obermann 


11 


52.6 


41-67 




Idiopathic 


5.5 ± 4.3 


11 


52.8 ± 11.6 


4/7 






Study 1 










blepharospasm 










2007 


Obermann 


9 


55.9 


43-63 


2/7 


Idiopathic cervical 


10 ± 6.8 


9 


57.5 ± 8.2 


2/7 






Study 2 








dystonia 










2006 


Rocca 


16 


42.7 


28-58 


1/15 


Aura 


24.8 


15 


38.6 24-50 


2/13 




2006 


Etgen 


16 


67.4 ± 4.3 




4/12 


Primary 


6.5 ± 4.9 


16 


65.3 ± 4.9 


4/12 
















blepharospasm 












2005 


Etgen Study 1 


28 


53.3 ± 8 




7/21 


Idiopathic restless 


19 ± 14 


28 


52 ± 8.3 


7/21 
















legs syndrome 












2005 


ttgen Study 2 


23 


59.3 ± 10 




6/17 


Idiopathic restless 


11 4 ± 11 


23 


59 ± 10.2 


D/ 1 / 
















legs syndrome 












2005 


Schimdt- 


40 


50.85 


16-70 


12/28 


Chronic tension 














Wilcke 








type headache; 
























medication- 
























overuse 
























headache 












2004 


Okada 


16 


34 


24-46 


10/6 


Chronic fatigue 


10 ± 244 m 


49 


34.4 21-47 


27/22 
















syndrome (CFS) 












2004 


Garraux 


36 


53 ± 9.7 




21/15 


Focal hand dystonia 


13 ± 7 


36 


52 ± 9.6 


21/15 














(FHD) 










1999 


May 


25 


47 


25-74 


23/2 


Idiopathic 




29 


33 20-55 


29/0 
















headache 











2.5. Distance matrix 

We then applied representational similarity techniques 
(Kriegeskorte et al, 2008) to the results obtained with network 
decomposition in order to investigate the relationship between 
networks. 



To obtain the visual representation, we first transformed the mod- 
eled activation maps of each experiment into a series of vectors con- 
taining the number of voxels altered in each network. Next, we con- 
structed a representational similarity matrix by calculating the cor- 
relation between each vector (r). The distance matrix (or represen- 
tational dissimilarity matrix) was derived as 1 - r. Lastly, we built a 
dendrogram of this matrix using a hierarchical clustering algorithm. 
This clustering approach had the advantage of allowing the visualiza- 
tion of hierarchies within the data. 
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2.6. Model-based clustering 

Studies about gray matter modifications often report similar 
groups of targeted areas. For instance area a is usually shown to be 
structurally altered together with area b. This raises the question of 
whether such modifications are observed together by chance or con- 
stitute a cluster. In this latter case, it may be that the alteration of one 
area is strictly related to the alteration of another. To investigate this 
possibility, we applied model-based clustering. 

The model-based clustering methodology is based on probability 
models, such as the Mixture Gaussian Model (Neumann et al, 2008). 
The idea behind this approach is that since a signal can be represented 
as a sum of sinusoids, each having its own amplitude, a probability 
distribution can be explained as a sum of Gaussians, each having 
its own weight. In mathematical terms, Mixture Gaussian models 
assume that the data is generated from a set of K simple Gaussian 
distributions like: 

k 

where 7T/< are the weights that a sample is drawn from, k is the 
mixture component and /x^and E k are the means and covariance 
matrices of a multivariate normal distribution N. Given a sample 
D = {X1.X2.X3. . . x n } the goal is to estimate the parameters of the 
mixture (7T/ 0 /z /<,£/<) for k = 1 . . . K. In other words, the idea is to find 
those parameters that best describe the original distribution, in terms 
of weights, mean and variance. In order to have a measure of how 
well the reconstructed distribution approximated the original one, 
we used the maximum likelihood estimation (MLE) approach that 
uses the logarithm of the likelihood function (L) for estimating pa- 
rameters: 

In P (D \7t, ii, E) = J2 ln J2 P* ( Zi = k ) Pn,E ( Xi 1 Zi = V- 

i k 

The MLE was used to identify the number of correct clusters us- 
ing as measure the Akaike information criterion (AIC) that is defined 
mathematically as: 

AIC = 2k-21n(L) 

where k is the number of parameters of the model (i.e. the number 
of clusters). The idea is to test a different number of clusters, (in 
our case from 1 to 20) and choose the number that minimizes the 
Akaike information criterion. The Akaike information criterion (AIC) 
constitutes a measure of the quality of the model and allows model 
selection. It is used to measure the goodness of fit of the model in 
relation to its complexity. 

2.7. Resting-state connectivity maps 

We observed two clusters of gray matter alterations in the anterior 
insula: one increase and one decrease (see also the Results section). 
Crucially, these clusters were located in two distinct anatomical po- 
sitions. We decided to investigate whether these two clusters were 
characterized by different connectivity patterns, using seed-based 
resting-state connectivity analysis. After acquiring resting-state data 
from twenty volunteers, we centered two ROIs in the two insular 
blobs and calculated their connectivity maps (see next paragraph for 
details). 

2.8. Resting-state data: subjects 

Twenty healthy right-handed volunteers (10 females, mean age 
32.6 ± 11.2) free of neurological or psychiatric disorders, not taking 
any medication known to alter brain activity, and who had no history 
of drug or alcohol abuse participated in this experiment. Handedness 
was ascertained with the Edinburgh Inventory (Oldfield, 1971). The 



written informed consent of each participant was obtained, in accor- 
dance with the Declaration of Helsinki. The study was approved by 
the institutional committee on ethical use of human subjects at the 
University of Turin. 

2.9. Task and image acquisition 

Images were acquired during a resting-state scan on a 1.5 Tesla 
Ingenia™ scanner (Philips Medical Systems). Functional T2* 
weighted images were acquired using echo planar (EPI) sequences, 
with a repetition time (TR) of 2000 ms, an echo time (TE) of 50 ms, 
and a 90° flip angle. The acquisition matrix was 64 x 64, with a 
200 mm field of view (FoV). A total of 200 volumes were acquired, 
with each volume consisting of 19 axial slices; slice thickness was 
4.5 mm with a 0.5 mm gap; in-plane resolution was 3.1 mm. Two 
scans were added at the beginning of functional scanning to reach 
steady-state magnetization before acquiring the experimental data. 
A set of three-dimensional high-resolution Ti -weighted structural 
images was acquired, using a Fast Field Echo (FFE) sequence, with 
a 25 ms TR, an ultra-short TE, and a 30° flip angle. The acquisition 
matrix was 256 x 256, and the FoV was 256 mm. The set consisted 
of 160 contiguous sagittal images covering the whole brain. 

2.10. Data analysis 

The data-sets were pre-processed and analyzed using the Brain- 
Voyager QX software (Brain Innovation, Maastricht, The Netherlands). 
Functional images were pre-processed as follows to reduce artifacts 
(Miezin et al, 2000): (i) slice scan time correction was performed 
using a sine interpolation algorithm; (ii) 3D motion correction was 
applied: using a trilinear interpolation algorithm, all volumes were 
spatially aligned to the first volume by rigid body transformations 
and the roto-translation information was saved for subsequent elab- 
orations; and (iii) spatial smoothing was performed using a Gaussian 
kernel of 8 mm FWHM. Several nuisance covariates were regressed 
out from the time courses to control for the effects of physiological 
processes (Bandettini and Bullmore, 2008; Birn et al., 2008; Napadow 
et al, 2008) and motion. Specifically, we included 2 additional covari- 
ates from white matter (WM) and cerebrospinal fluid (CSF), as well as 
6 motion parameters. The time courses were then temporally filtered 
in order to keep frequencies between 0.008 and 0.08 Hz only, and nor- 
malized. Subsequently, the datasets of each subject were transformed 
into Talairach space (Talairach and Tournoux, 1988). Two spherical 
ROIs with a diameter of 8 voxels were placed in the center of mass 
of the two insular blobs derived from the results of the model-based 
clusterization. The coordinates of the ROIs of increased gray matter 
were 31,11,4; those of the ROIs of decreases in gray matter were -30, 
10, -5. For each subject the time courses of the two ROIs were ex- 
tracted and normalized. A multi-subject random effect General Linear 
Model was performed using each subject's time courses as predictors. 
All maps were thresholded at p < 0.05 and corrected for multiple 
comparisons using the false discovery rate (FDR). 

3. Results 

3.1. Anatomical likelihood estimation 

The ALE results showed that chronic pain is associated with a com- 
mon core set of gray matter decreases in the bilateral medial frontal 
gyri, bilateral superior frontal gyri, right pre- and post-central gyri 
(including the primary somatosensory and primary motor cortex), 
bilateral insula (anterior), right cingulate cortex (dorsal posterior cin- 
gulate cortex), basal ganglia, thalamus and periaqueductal gray. In- 
creased gray matter was found in the bilateral post-central gyrus, left 
inferior parietal lobule, right pre-central gyrus (primary motor cor- 
tex), right post-central gyrus, in the dorsal prefrontal areas, in the 
caudate, thalamus, cerebellum, and pons (see Fig. 1 and also Tables 
2-3). 
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Table 2 

Gray matter decreases: activation likelihood estimation results. 



Cluster # 


Volume (mm 3 ) 


Weighted center (x, y, z) 




X 


y 


z 


Label 


1 


1368 


58.5 


-9.51 


12.97 


60 
58 


-8 
-12 


10 

22 


Right precentral 
gyrus BA 43 
Right postcentral 
gyrus BA 3 


2 


808 


11.53 


-23.63 


11.81 


12 


-24 


12 


Right pulvinar 


3 


672 


22 77 


14 62 


-5.36 


20 
28 


12 
20 


-6 
-6 


Right putamen 
Right insula BA 47 


4 


608 


-4.72 


36.35 


26.03 


-6 


36 


28 


Left medial frontal 
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Table 3 

Gray matter increases: activation likelihood estimation results. 
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Fig. 1. Gray matter anatomical likelihood estimation (ALE) results. This image sum- 
marizes the results of all the papers involved in this meta-analysis. Colors from red to 
yellow show gray matter increases, colors from blue to green show gray matter de- 
creases (two-dimensional ALE maps were computed at a false discovery rate corrected 
threshold of p < 0.05, with a minimum cluster size of K > 100 mm 3 and visualized 
using MRIcron). Talairach coordinates of clusters showed in this images are reported 
in Tables 2 and 3. 




Fig. 2. Results of the jackknife analysis on the ALE GM decrease results. A high prob- 
ability (approaching 100%) means that a certain area is led by the majority of the 
experiments included. By contrast, a lower probability means that a certain area is led 
by few experiments thus indicating the need to exercise caution in interpreting this 
point. Talairach coordinates of clusters showed in this images are reported in Table 2. 




Fig. 3. Results of the jackknife analysis on the ALE GM increase results. A high prob- 
ability (approaching 100%) means that a certain area is led by the majority of the 
experiments included. By contrast, a lower probability means that a certain area is led 
by few experiments thus indicating the need to exercise caution in interpreting this 
point. Talairach coordinates of clusters showed in this images are reported in Table 3. 



3.2. Jackknife analysis 

The jackknife analysis (Tukey, 1958) (Figs. 2 and 3) showed all the 
areas found in the ALE maps to have very high reliability. Areas with 
lower (although still very high) reliability were found in the thalamic 
and ventral prefrontal cortices. Overall the reliability of gray matter 
increases was 20% lower than that of decreases, probably due to the 
paucity of papers reporting gray matter increases. 

3.3. Network decomposition 

With the advent of models allowing the investigation of functional 
integration rather than mere functional segregation (see Friston, 2009 
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Fig. 4. Involvement of the large-scale brain networks in chronic pain-related gray 
matter alterations/The upper panel shows the global involvement of each large-scale 
brain network in the pain-related gray matter alterations (the results from all the 
pathologies are here placed together, the values are expressed in mm 3 of altered brain 
volume). All the networks except auditory, VI , V2 and V3 are significantly different from 
zero. The lower panel shows a cumulative graph depicting the different involvement of 
the large-scale brain networks in different chronic pain pathologies (the results from 
all the pathologies are split here, the values are expressed in mm 3 of altered brain 
volume). To obtain this graph the data were centered (the mean was subtracted). So 
higher values represent a greater deviation from the mean. 



and also Farmer et al, 2012), researchers have begun to use a 'network 
approach' and study the co-variation of brain activity in response to 
noxious stimuli in order to unveil the functional significance of brain 
responses to those stimuli (Farmer et al, 2012). This approach is able 
to provide new insights in the study of the complex mechanisms 
underpinning the emergence of pain (Davis, 2011). To investigate 
which of the large-scale brain networks show GM alterations, we 
compared the number of increased/reduced voxels in each of the 
brain networks described by Biswal et al. (2010). 

As shown in the upper left panel of Fig. 4, each type of chronic pain 
shows a different involvement of each large-scale network. To obtain 
this graph the data were centered (the mean value was subtracted); 
as a result, lower values represent a small deviation from the mean. 
The maximal variability was expressed by the thalamic-basal ganglia 
network (Th-Ba), followed by the DMN and premotor and somatosen- 
sory networks. Trigeminal pain showed a reduced involvement of the 
Th-Ba and DMN networks whereas complex regional pain syndrome, 
blepharospasm, chronic fatigue syndrome and back pain were charac- 
terized by reduced involvement of the somatosensory and premotor 
networks. Interestingly, the salience and attentional networks were 
damaged in a very similar way by different chronic pain pathologies. 

The lower panel of Fig. 4 shows the mean number of voxels 
involved for each network. The DMN and Th-Ba, attentional and 
salience networks are the areas altered most. 
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Fig. 5. Representation of dissimilarity between each large-scale brain network. The 
upper panel shows the (dis)similarity between each network. Networks that are in- 
volved in the same kind of pathologies are in a closer position. Networks that are 
involved in different pathologies are placed apart.The dendrogram in the lower part 
shows the results of a hierarchical clusterization of the large-scale brain networks on 
the basis of the different involvement of each network. Networks that show a similar 
involvement in the same kind of pathologies are clustered together. 




3.4. Distance matrix 

The (dis)similarity representation analysis (Fig. 5) showed that 
with the exception of the thalamus-basal ganglia network, the other 
networks could be clustered into three groups showing a similar in- 
volvement in chronic pain pathologies: a first group was composed 
of the DMN and cerebellar and motor networks; a second group of 
the salience and attentional networks; and a third group of the OFC, 
premotor, sensorimotor, auditory and visual networks. 




Fig. 7. Resting-state seed voxel functional connectivity of the two anterior insular gray 
matter alterations: decrease-related and increase-related. The lower left panel shows 
the areas functionally connected to the "increase seed". The lower right panel shows 
the areas functionally connected to the "decrease seed". 



3.6. Resting-state connectivity maps 



3.5. Model-based clustering 

The optimal number of clusters was calculated by minimizing the 
Akaike information criterion. This method showed an optimal cluster 
number of three for gray matter decreases and two for gray matter in- 
creases. The left panel of Fig. 6 shows the three clusters into which the 
gray matter decreases were decomposed. The first cluster was found 
to involve the fronto-parietal and medial wall areas, the second in- 
cluded the operculo-insular, cingulate, posterior thalamic and medial 
prefrontal areas, and the third comprised the temporal and pontine 
areas. The right panel of Fig. 6 shows the three clusters into which the 
gray matter increases were decomposed. The first included the so- 
matomotor and somatosensory, premotor and parietal areas and the 
second the operculo-insular, basal ganglia, pontine and dorso-ventral 
prefrontal areas. Interestingly both decreases and increases were ob- 
served in the anterior insula, but in different portions of the region 
(see Fig. 7). 



We performed a resting-state connectivity experiment to explore 
whether the two anterior insular clusters derived from the decompo- 
sition of gray matter increases and decreases belong to different net- 
works. As shown by Fig. 7, the two insular clusters presented a very 
different functional connectivity pattern. The increase cluster had a 
strong connectivity to the saliency detection network with anticor- 
relations with a series of medial, dorsal prefrontal and parietal areas 
of the default mode network (DMN). The decrease cluster showed 
a connectivity with areas anticorrelated with the increase cluster. 
These areas were mainly constituted by the anterior part of the DMN. 
The decrease cluster showed anticorrelation with the salience detec- 
tion network. These results suggested that the functional connectivity 
profile of the two ROIs is coherent with the gray matter modifications. 
The opposite alteration (increase or decrease) was mirrored in the op- 
posite functional connectivity of these two ROIs: while the "increase" 
ROI showed a correlation with the salience network and an anticorre- 
lation with the DMN, the "decrease" ROI presented an anticorrelation 
with the salience network and a correlation with the DMN. 
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4. Discussion 

This study was designed to: i) verify the presence of a core set of 
brain areas commonly modified by chronic pain; ii) investigate the 
involvement of these areas in a large-scale network perspective; iii) 
study the relationship between altered networks and; iv) find out 
whether chronic pain targets clusters of areas. Our results show that: 
i) gray matter alterations in single areas should be better conceived 
in a framework of large-scale brain networks and ii) large-scale brain 
networks present both pathology-unspecific and pathology-specific 
involvements. 

4. 1 . Network approach for the study of gray matter modifications 
induced by chronic pain 

It has been demonstrated that the presence of long-lasting ongoing 
pain can modify the structure of the brain, inducing local morpholog- 
ical changes in the brain parenchyma (Farmer et al, 2012;May, 2008, 
201 1 ; Woolf and Salter, 2000). Indeed, in the last decade, the applica- 
tion of neuroimaging techniques, such as voxel-based morphometry 
(VBM), has provided considerable insight into structural brain reor- 
ganization in subjects suffering from chronic pain syndromes. The 
results of previous VBM studies have converged in concluding that 
chronic pain induces gray matter structural changes, often related to 
the duration of pain (Apkarian et al, 2009). However, these studies 
have not always obtained homogeneous results as to the areas af- 
fected by chronic pain and the amount of overall reorganization. In 
the present study we characterized gray matter modifications under a 
novel perspective which took into account networks rather than sin- 
gle areas. Our choice was motivated by recent theoretical frameworks 
which underline the importance of large-scale networks in the elabo- 
ration of noxious stimuli (Cauda et al, 2013; Mayhew et al, 2013), in 
the transition from acute to chronic pain (Baliki et al, 2012; Hashmi et 
al, 2013) and in the maintenance of chronic pain (Farmer et al, 201 2). 
Indeed, it is now widely accepted that brain regions work in synergy 
and that their correlated or anticorrelated activity can be grouped into 
several large-scale networks (Fox et al, 2005; Mesmoudi et al, 2013). 
Besides, single areas can belong to more than one network and be 
functionally linked to different regions depending on the task (Cauda 
et al, 2012; Torta et al, 2013). This network approach may help to 
explain why previous studies on gray matter alterations in chronic 
pain have reported inconsistent results at the level of single areas. 

4.2. Shared and specific alterations in large-scale brain networks 
induced by chronic pain 

Our results revealed that chronic pain can induce both common 
and pathology-specific changes in large-scale networks. Gray mat- 
ter alterations were almost always found to be present, regardless 
of the chronic pain syndrome, in the DMN and the thalamus-basal 
ganglia circuit and attention networks. In contrast, sensory networks 
were seldom and variably targeted by chronic pain. Trigeminal, mus- 
culoskeletal, phantom limb pain and irritable bowel syndrome were 
characterized by a minor involvement of the DMN and basal ganglia 
networks, but by greater alterations of the somatosensory and pre- 
motor networks. Conversely, blepharospasm, complex regional pain 
syndrome and chronic low back pain were characterized by a greater 
involvement of the DMN and basal ganglia networks, with a more 
marginal contribution of the somatosensory and premotor networks. 
Baliki et al. (201 1 ) have recently investigated gray matter alterations 
in chronic back pain, osteoarthritis and complex regional pain syn- 
drome considering gray matter alterations under an area and network 
perspective. The authors compared the three groups and a healthy 
control group in terms of total neocortical volume reduction and ev- 
idenced that only the chronic back pain group differed significantly 



from the others. In addition, by performing a voxel based morphom- 
etry analysis, the authors observed common and specific reductions 
in gray matter for the three pain conditions. All groups showed de- 
creases in the insula. The chronic back pain group had reductions in 
somatosensory areas, pre- and post-central gyri, the hippocampus 
and temporal lobes; the osteoarthritis group in the cingulate cortex, 
hippocampus and inferior temporal cortex; and the complex regional 
pain syndrome group in the orbitofrontal cortex. Moreover the au- 
thors parceled the cortex into 82 ROIs and observed that secondary so- 
matosensory cortices, bilateral insulae and dorsal and orbital frontal 
regions, besides the hippocampus were those that better differenti- 
ated between groups. In order to study gray matter alterations in an 
inter-relation between areas, the authors correlated the gray matter 
decreases between each pair of ROI: they observed specific alterations 
in pain patients as compared to healthy subjects. Specifically, whereas 
healthy controls showed greater correlations between neighboring 
voxels (i.e. anatomically contiguous areas were more likely to have 
similar gray matter density), pain patients presented strengthened 
long distance correlations. These alterations were specific enough to 
be able to classify the group correctly. 

Ours is however the first report that, by covering a wide spec- 
trum of chronic pain syndromes, is able to specify that chronic pain 
targets large-scale brain networks. This is in line with the network de- 
generation hypothesis characterized in degenerative syndromes that 
posits that atrophy (here conceived as gray matter decreases) reca- 
pitulates healthy brain functional networks (Seeley et al, 2009). In 
this view, the pathology targets specific large-scale networks, which, 
in the healthy brain, feature the perceptual or cognitive function that 
is altered in the disease. 

Notably, our findings show that attention circuits (the salience, 
ventral and dorsal attentional networks) are similarly affected in all 
pain syndromes and their involvement is not subject to change in re- 
lation to the pathology under examination. This finding can substan- 
tiate the idea that attention networks are always targeted by chronic 
pain, regardless of its etiology and regardless of the attention network 
(see Cauda et al, 2010). Indeed, we found gray matter decreases to 
affect both the ventral exogenous attention network, usually linked 
to bottom-up processing, and the dorsal endogenous attention net- 
work, linked to top-down processing. We also observed that the DMN 
was targeted by all chronic pain pathologies. The DMN is thought to 
be implicated in self-referential cognitive processing, self-awareness 
and self-monitoring. It is plausible to hypothesize that a long-lasting 
and intruding disorder like chronic pain alters interoceptive functions 
per se, regardless of the etiology of the disease. 

The thalamic-basal ganglia networks were also always targeted. 
Thalamic-cortical dysrhythmia has often been proposed as a mech- 
anism that underlies chronic pain states (Llinas et al, 1999). The ob- 
servation of increased levels of power spectrum in various pain states 
measured by EEG (Saab, 201 2 ; Sarnthein et al, 2006; Stern et al, 2006) 
has led to the suggestion that altered brain oscillation may constitute 
a characteristic of dysfunctional thalamo-cortical loops. Recently, it 
has been shown that, in chronic neuropathic pain, the thalamus is 
characterized by significant volume loss and reduced inhibitory neu- 
rotransmitter content (Henderson et al, 2013). 

4.3. Investigating how chronic pain affects large-scale brain networks 

By applying representational similarity techniques to the re- 
sults of the network decomposition we observed a panellation into 
three clusters (Fig. 5 lower panel), where the first cluster was com- 
posed of the DMN, cerebellum and motor areas, the second of the 
salience/attention networks and the third of the sensory networks. 
The thalamic-basal ganglia network remained isolated. 

These results substantiate our previous conclusion: that the atten- 
tion and sensory networks are differentially affected by chronic pain 
pathologies: the attentional networks show pathology-unspecific 
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changes, whereas in the sensory networks strong damage variabil- 
ity is observed for different chronic pain categories. 

To obtain true data-driven confirmation as to whether chronic pain 
specifically targets some large-scale brain networks we performed a 
model-based clusterization of our data-set. In this analysis we inves- 
tigated increases and decreased in gray matter separately. 

Our analyses confirmed the network-based results, indeed this 
new data-driven analysis highlighted three clusters of decreases and 
two of increases. Increases were observed not only in the somatosen- 
sory and somatomotor areas, as previously reported (Schweinhardt 
et al, 2008), but also in the operculo-insular cortex. Decreases were 
observed in the fronto-parietal regions, temporal and pontine regions 
and, again, in the anterior insula. The fact that decreases and increases 
tap on differential networks is important. As already suggested by 
previous studies on chronic pain and now supported by data on ex- 
perimental pain, it is highly likely that differential networks respond 
to noxious and painful stimuli in a different manner and with dif- 
ferent temporal envelopes (Cauda et al, 2013; Mayhew et al, 2013; 
Moulton et al, 2012). These networks may thus undergo selective 
alterations by the pathology. 

4.4. Role of the insular cortex in gray matter alterations 

Our findings indicate that both increases and decreases in gray 
matter can be found in the insula. Crucially, increases and decreases 
did not occur in the same region of the anterior insula: increases 
were observed in the dorsal anterior insula, whereas decreases were 
detected in the ventral anterior insula. These two areas, as revealed 
by resting-state functional connectivity, are correlated with different 
networks. The ventral anterior insula (gray matter increases) was cor- 
related with the salience network and anticorrelated with the DMN. 
Conversely, the dorsal anterior insula was correlated with anterior 
parts of the DMN and anticorrelated with the salience network. These 
data offer another possible explanation for the increase /decrease pat- 
tern that has so far been explained in terms of duration of the pathol- 
ogy (Apkarian et al, 2004). Our results suggest that increases and 
decreases occur in two anticorrelated networks: increased activity 
in one network leads to decreased activity in the other. This view in 
which different areas of the brain that are related to different and 
somewhat opposite functions exert an opposite (anticorrelated) con- 
nectivity pattern and are attacked in an opposite way by pathologies 
opens novel perspectives (Mesmoudi et al, 2013). 

A possible explanation for the pattern of decrease/increase in the 
insula may be that the salience network and sensory regions are over- 
excited by greater weights attributed to external stimuli such as pain. 
This functional over-recruitment may promote a structural reorgani- 
zation of the areas in the network (such as the somatosensory areas 
and ventral anterior insula) and induce an increase in gray matter. 
The excessive activation of the salience network may, in turn, cause 
de-activation of the anticorrelated circuit (DMN), which, as a conse- 
quence, undergoes a reduction in volume. This explanation implies a 
tight link between functional and structural modifications, which has 
yet to be fully characterized in chronic pain. However, at least two 
lines of reasoning support this view. First, it is difficult to conceive 
that functional and structural alterations are not related to one an- 
other since functional and structural networks are (Kelly et al, 2012; 
Pernice et al, 201 1 ; Zhang et al, 201 1 ). Numerous fMRI studies have 
already demonstrated that chronic pain alters the normal functioning 
of brain networks (Balenzuela et al, 2010;Cauda et al, 2009a 2010). 
For instance, it has been shown that chronic pain disrupts thalamo- 
cortical connections (Cauda et al, 2010; Jensen et al, 2012), default 
mode network dynamics (Baliki et al, 2008 ; Napadow et al, 2012) and 
the functional connectivity of the dorsal and ventral attentional net- 
works (Cauda et al, 2010; Napadow et al, 2010). Second, it has been 
shown that synchronous neuronal firing promotes synaptogenesis 
(Bi and Poo, 1999; Katz and Shatz, 1996) and that intrinsic functional 



connectivity and structural covariance converge (Seeley et al, 2009). 

Many areas of decreases in gray matter that we found in our data 
belong to the descending inhibitory circuits. Therefore, another possi- 
ble conclusion is that the unbalance between anticorrelated networks 
in favor of increased salience towards external stimuli hampers the 
physiological inhibition of potentially dangerous stimuli. This mech- 
anism may be another candidate for the maintenance of chronic pain 
states. This explanation is supported by theoretical models which 
suggest a less efficient top-down modulation of pain in chronic pain 
patients. In this sense, chronic pain patients may be more prone to 
attentional processing of pain-related information (see Van Damme 
et al, 2010 for a critical review on attentional-related aspects of pain 
processing) and might have difficulties in disengaging attention from 
pain (Van Damme et al, 2004). 

4.5. Limitations and future perspectives 

As our analyses relied on data provided by previous studies, we 
were not able to control for factors that could have intervened in 
altering gray matter and that had not been tested in the original re- 
search. This is important, as chronic pain patients frequently have 
several comorbid conditions, for example anxiety and mood disor- 
ders, altered life-styles and drug use that might directly contribute to 
morphological changes. These aspects have not always been consid- 
ered in previous VBM studies. In particular, the presence of chronic 
anxiety, depressed mood and chronic stress leads to activation of the 
hypothalamic-pituitary-adrenal (HPA) axis that, by increasing stress 
hormones (Cortisol in humans), induces neuronal loss in several brain 
areas. This loss is related to several factors such as an increase in pro- 
inflammatory cytokines, a reduction in neurotrophic factors and a 
rise in neurotoxic glutamate (Swaab et al, 2005; Tata and Anderson, 
2010). Chronic pain shares several of these mechanisms (increase 
in cytokine and glutamate release) that may induce a vicious cycle 
in perpetuating the comorbidity between pain and emotional disor- 
ders and between the latter and neuronal loss (Torta and Munari, 
2010). It has been observed that, in fibromyalgia patients, the volume 
decrease in regions that process emotional states is related to the pa- 
tient's affective state and not to the presence of pain per se (Hsu et 
al, 2009). What is more, the chronic use of opioids produces neuro- 
plastic changes in animals and also in humans, as confirmed by gray 
matter dose-related volumetric changes (Younger et al, 201 1 ). In line 
with this, the use of antidepressants in chronic pain can offer pro- 
tection against neuronal suffering, by preventing the enhancement of 
glutamate release induced by stress (Musazzi et al, 2011). 

In conclusion, our results showed that chronic pain selectively al- 
ters large-scale brain networks. Some of these networks are altered 
in a very similar way but others are altered differentially by chronic 
pain pathologies. This adds to a picture in which some common areas 
represent a "core group" of regions altered by almost all the chronic 
pain pathologies and some other areas that are differentially dam- 
aged and that may represent the specific damage of each pathology. 
Indeed as already proposed by Baliki et al. (201 1 ), each syndrome also 
has a specific and discriminable pattern of alterations. Future studies 
including a larger sample will help to establish the validity of such 
suggestions. 
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